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The death of FRII radio sources and their connection with radio 
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ABSTRACT 

Radio relic sources in galaxy clusters are often described as the remnants of powerful radio 
galaxies. Here we develop a model for the evolution of such relics after the jets cease to 
supply energy to the lobes. This includes the treatment of a relic overpressured with respect to 
its gaseous surroundings even after the jets switch off. We also determine the radio emission 
of relics for a large variety of assumptions. We take into account the evolution of the strength 
of the magnetic field during the phase of relativistic particle injection into the lobes. The 
resulting spectra show mild steepening at around 1 GHz but avoid any exponential spectral 
cut-offs. The model calculations are used to fit the observed spectra of seven radio relics. The 
quality of the fits is excellent for all models discussed. Unfortunately, this implies that it is 
virtually impossible to determine any of the important source parameters from the observed 
radio emission alone. 

Key words: galaxies: active - galaxies: jets - intergalactic medium - radio continuum: galax- 
ies 



1 INTRODUCTION 

Diffuse radio emission without an apparent host galaxy is found in 
an increasingjiumber^ofcl^^ et al., 



1999; ISiovannini & Feretti, 200Cl|;|Kempner & Sarazin, 200 1|; G ov- 
oni et al., 2001). These objects can roughly be divided into radio ha- 
los close to the cluster centre which appear comparatively smooth 
in radio images and radio relics with a far more distorted and knotty 
appearance. The latter are usually fo und further away from the cen- 
tral regions ( Slee & Reynolds, 1984 ). Here we concentrate on radio 



relic sources. 

A large number of possible explanations for radio relics has 
been proposed. All of these aim at explaining the production of a 
population of relativistic electrons responsible for the observed ra- 
dio synchrotron emission. The suggested mechanisms range from 
the turbulent wakes of galaxies moving in the cluster's gravitational 
potential to cluster mergers. In the present paper we concentrate 
on the picture of relics as remnants of once powerful radio galax- 
ies or radio-loud quasars. This scenario is based on the idea that 
the relics are the lobes of former radio galaxies, the jets of which 
ceased to supply themwithenerg^soi^^ 

& Gubanov, 1994; [Slee et al., 200 1[ ). The relativistic electrons, the 
emission of which we observe today, were accelerated at the strong 
shocks at the ends of the active jets. The electrons lose their energy 
due to synchrotron losses, inverse Compton scattering of cosmic 
microwave background photons and, possibly, further adiabatic ex- 
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pansion of the relics themselves. As the relativistic electrons are not 
replenished by the active jets, the radiative energy losses introduce 
a spectral cut-off moving towards low frequencies in time. This can 
potentially explain the observed steep radio spectra of relics. 

Note however that not all relics may be consistent with the 
picture of the passively fading remnant of a radio galaxy. Some 
relics are extended and show a rather smooth structure resembling 
in some ways the appearance of radio haloes, while others are com- 
posed of a number of individual knots. The variety of morphologies 
different from the appearance of active radio galaxies may very 
well indicate that our interpretation of relics as the remnants of ra- 
dio galaxies is incorrect in many cases. However, it is not clear what 
the fate of the material forming the large scale radio structure of a 
formerly active radio galaxy is. Buoyant rise in the gravitational 
potential well of the cluster will certainly play a role (e.g. Chura- 
zov et al., 2001), but the movement of the parent galaxy through 
the cluster gas and cluster mergers may fragment the remnant as 
well. It is therefore not surprising that the morphology of relics 
is seldomly reminiscent of active radio galaxies. Even if all relics 
were indeed, as we assume here, the remnants of radio galaxies, the 
evolution of the remnant may imply significant departures from the 
simple pressure evolution envisaged in the model introduced here. 

Another problem is the possibility of re-acceleration of the 
relati vistic electrons even after the jets have switched off. S lee et al. 
(2001) present radio observations of the relics in A13, A85, A133 
and A4038. The substantial substructure in these sources may well 
be the signature of some highly dynamic processes being at work 
which could lead to further acceleration of relativistic particles. In 
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this paper we will neglect any re-acceleration processes which may 
take place in relics. 

The question then is whether or not a passively evolving popu- 
lation of relativistic electrons produces a radio spectrum consistent 
with observations. The main problem is that the radio spectra of 
relics are shallow at around 100 MHz and considerably steeper in 
the GHz; but they do not show exponential cut-offs. The spectra 
are therefore not consistent with a single population of electrons 
accelerated all at the same time and passively losing their energy. 
Taking this ext ended injection period into account, K omissarov & 
Gubanov (1994) showed that two breaks appear in the spectrum. 
The lower break is located at the frequency cut-off of the 'oldest' 
electrons while the higher break is determined by the 'youngest' 
electrons. The comparatively mild steepening of the spectra be- 
tween the breaks is roughly consistent wit h the observations. Nev- 
ertheless, Komissarov & Gubanov (1994) found that for effective 



pitch-angle scattering of the relativistic electrons the steepening of 
the spectrum was still too strong. Therefore they had to suppress 
pitch-angle scattering to further flatten their model spectra. 

The exponential cut-off in the spectrum is also avoided in the 
case of inhomogeneous magnetic field strengths inside the relic. 
The relativistic electrons spend most of their time in the low-field 
regions and only occasionally diffuse into the high-field regions. 
Thus their energy losses are reduced and by careful adjustment of 
the efficiency of the diffusion process an emission spectrum with 



mild steepening at high frequencies c an be achieved ( Fribble, 1994 



Eilek et al., 1997). Slee et al. (2001) fitted a model based on this 



idea to the spectra of relic sources a nd found good agreement. 

Komissarov & Gubanov (1994) assumed that the relics or their 
progenitors were in pressure equilibrium with their surroundings 
even at the time of the injection of relativistic electrons. In this pa- 
per we expand this model by taking into account the evolution of 
the lobes during the injecti on of relativistic particle s. If the lobes of 
radio galaxies of type FRII ( Fanaroff & Riley, 1974 ) are the progen- 



itors of the radio relics, then during the time the jets are still active, 
the lobes are expanding into the surrounding material. This implies 
that the pressure inside the lobes and therefore also the strength of 
the magnetic field are not constant during the particle injection. In 
fact, the strength of the magnetic field is decreasing and thus rela- 
tivistic electrons injected at later times can survive for longer. We 
show that this effect results in radio spectra of the resulting relic 
sources without exponential cut-offs and a spectral slope compara- 
ble to that in observed relics. 

At the time the jets switch off, the lobes may still be over- 
pressured with respect to the ambient gas and therefore continue to 
expand. We derive expressions for the temporal behaviour of the 
volume and the pressure of such 'coasting' relics. We take such a 
coasting phase into account in our calculation of model radio spec- 
tra from relic sources. Finally, we show that our models fit the ob- 
served spectra very well. Unfortunately, we find that none of the 
important source parameters, like the source age, or the physical 
processes taking place in the relics can be determined from the ra- 
dio spectra alone. 

In Section |^ we derive the temporal behaviour of a coasting 
relic source after the jets have switched off. We use the results 
found to calculate the radio spectra of relic sources in Section |3| 
Here we also show how the variation of the main model parameters 
influences the results. In Section ^ we fit the observed spectra of 
relic sources with our model and show that very little useful infor- 
mation about the conditions in the sources and their environments 
can be gained from the models. We summarise our results in Sec- 
tion Isl 



2 THE SPHERICAL RADIO SOURCE 

After the jets of a powerful radio sources switch off, the cocoon 
inflated during the active phase may still be overpressured with re- 
spect to the external medium. It will then enter the 'coasting phase' 
while continuing to expand behind a strong bow shock. This situa- 
tion is somewhat reminiscent of the early evolution of a supernova 
remnant or the explosion of a nuclear weapon. The dynamics of a 
strong explosion in a uniform atmosphere are described by the sim- 
ilarity solutioi^^of^ylor^an^^Sedo^^(e^gJSedw & 
Lifshitz, 1987). It is straightforward to generalise this solution for 
the case of a power-law density distribution in the external medium 
with p = po(r/ao)^P, where ay is the scale height or core radius 
of the distribution. We will assume this density distribution for all 
model calculations in this paper. The radius of the spherical bow 
shock then grows as i?s oc f2/(5-|3) pgspifg the fact that the co- 
coons of radio galaxies are usually not spherical, one might expect 
that their expansion during the coasting phase is governed by the 
same proportionality once the jet activity has ceased. However, in 
the following we will show that the presence of the contact discon- 
tinuity defining the cocoon boundary in radio galaxies requires a 
different dynamical behaviour. 

Goldshmidt & Kephaeli (1994) present a different solution to 
this problem based on their assumption that the dynamics of the ra- 
dio source expansion are not influenced by the mass of the external 
gas swept up by the bow shock. The dynamical behaviour of active 
radio sources strongly suggests that thei r cocoons ar e underdense 
with respect to the external medium (e.g. Falle, 1991). Therefore it 



is very unlikely that this assumption can be made for the problem 
considered here. 



2.1 Analytical solutions 

For simplicity we will assume that the bow shock and the cocoon 
are both spherical. In practice both are elongated along the jet axis 
but this different geometry results only in a change of the constants 
of proportionality. The dynamical behaviour, i.e. the time depen- 
dence of the physical parameters of the flow, will be the same. 
The cocoon volume is then given by Vc = 4/37ti?^ and that of the 
layer of shocked gas in between the bow shock and the cocoon is 
Vs = 4/3n(i?s — RI)- Furthermore, we also assume that the pres- 
sure, p, inside the cocoon and in the shocked layer surrounding 
it, is uniform. Again, this greatly simplifies the problem and will 
not change our results apart from the exact numerical values of 
the normalisations. During the active phase, while the jets are still 
switched on, conservation of energy for the gas inside Vc requires 



-urIp + ^hJcPRIRc 



(Yc-l)ej, 



(1) 



where Yc is the adiabatic index of the cocoon material and Qj is the 
energy transport rate of the jets. We will keep constant during 
the lifetime of the jets at the end of which it drops to zero instanta- 
neously. Similarly, for the shocked gas inside we find 



n(Rl-Rl^p + 47tYsP (R-IR, - R^R, 



27t(Ys-l)/f?Po 



Rl 



(2) 



where Ys is the adiabatic index of the shocked material. The ex- 
pression on the right of equation ^ is the rate at which kineti c 
energy flows into the bow shock (see [Kaiser & Alexander, 1997 ). 
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Finally, since we consider only strong shocks, the pressure within 
the source is balanced by the ram pressure of the external medium. 



(3) 



Equations through are identical to the system of equations 
studied by Reynolds & Begelman (199"^ for the case of intermittent 
jet sources. 

We are interested in steady-state similarity solutions for this 
set of equations with the three unknown functions i?c, and p 
being power laws of time. We find 



(Yc-1)(5-P)^ ,fei 
12jt(9Yc-4-p) 

3 \ l/(5-P) 



1/3 



p 



Po«c 



Poao 



-4-p 



3 \ l/(5-P) 



(4) 



l/(5-P) 



where 



(ye-l)(5-p)^(9Ys-4-p) 
67i(9Yc-4-P) (9Ys + l-2p)' 



(5) 



The exponen ts of the time vari able in these equations are e qual 
to those found by Falle (1991 ) and Kaiser & Alexander (1997 ) for 
sources with active jets. 

After the jets cease to supply energy to the cocoon, the source 
enters the coasting phase and we have to set = 0. We assume 
that the cocoon stays overpressured with respect to the surround- 
ing medium. In this case, equation (hi) describes the now adia- 
batic expansion of the cocoon, i.e. p °^ Rc ■ With this, equations 
(^ and (^) in principle again yield power-law similarity solutions 
with oc R^oc f2/(2-|3+3yc)^ j^^jg ^^^^ d^j^ different from the 
Taylor-Sedov solution. This difference is caused by the adiabatic 
behaviour of the cocoon driving the bow shock for which there is 
no analogon in the case of a strong explosion. The ratio of the ra- 
dius of the spherical cocoon and the bow shock is given by 



Rc 
R. 



1 



11/3 



2(Ys-Yc) 



(6) 



Obviously, this solution cannot be correct for Yc = Ys or Ys = 2Yc — 1 
which includes the case of a relativistic cocoon (Yc =4/3) sur- 
rounded by 'cold' material (Ys ~ 5/3). Nevertheless, we will see 



in section 2.2 that in the latter scenario the steady-state solution is 
an excellent approximation to the evolution of R^. 

For Yc = Ys all terms containing Rc in equation (^ cancel and 
a different steady-state solution is possible with 



^c 

Rs 



2(Tc + l) 
;Tc(7+3Tc-2P) 



-6(Tc- 



(7) 



p cx /7+3TC- 



This implies that the ratio Rc/Rs is no longer a constant. Despite 
this, the expansion of the bow shock and that of the cocoon are still 
intrinsically self-similar. 




t / yeors 

Figure 1. Comparison of the numerical results with the analytical solu- 
tions for Yc = Ys = 5/3. Upper panel, solid lines: Numerical results for 
Rc (lower line) and Rg (upper line); short dashed lines: analytical solution 
for active jets; long dashed lines: analytical solutions for 'coasting phase'. 
Lower panel: Analytical solutions divided by numerical solutions. Dashed 
line: Ratio for R^, solid line: Ratio for Rs. 



2.2 Numerical solutions 

All of the solutions derived in the previous section describe steady- 
state configurations. To assess how quickly a spherical radio source 
would change from one of these states to another when the jets are 
switched off, we also integrated the equations (pi) to (pl) numeri- 



1.7 X 10 



gem 



cally. For this we set Qj = 10 ergs s , po 
ao = 500 pc and P =1.5. The same values were used by Reynolds 
& Begelman (1997) and they imply that poflg = 10^ gcm"' '^. At 
a source lifetime of 10^ years, we 'switched off' the jets by set- 
ting Qj = 0. For the integration we used a standard Runge-Kutta 
algorithm. The initial conditions were derived from equations (^. 

Figure |l| shows the results of this numerical integration for 
the radius of the bow shock, R^, and that of the cocoon, R^ for the 
case Yc = Ys = 5/3. The analytical solutions for the still active jets 
are calculated from equations (^). For the coasting phase we used 
equations (Q). Here we are free to choose the normalisation of the 
analytical solutions as these are not constrained by the equations. 
We have chosen the normalisation such that the analytical solution 
agrees with the numerical results at 2 x 10^ years. The lower panel 
of Figure [l] shows the ratio of the analytical solutions and the nu- 
merical results. We have used the analytical solutions for the active 
phase for times less than 10^ years and those for the coasting phase 
for later times. The discrepancy between the two solutions is al- 
ways less than 20%. 

Figure ^ compares the numerical results for the case Yc =4/3 
and Ys = 5/3 with the steady-state similarity solutions. As ex- 
pected, the agreement is worse. Particularly the radius of the bow 
shock grows more quickly than predicted by the power-law solu- 
tion. However, the radius of the cocoon is still well described by 
the steady-state solution. The disagreement between the numeri- 
cal integration and the analytical power-law is comparable for this 
case to that for Yc = Ys- We therefore conclude that the dynamical 
behaviour of the cocoon of a spherical radio source adjusts quickly 
to the cessation of the energy supply by the jets. 
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Figure 2. Same as figure ^but witli = 4/3 and -/j = 5/3 

2.3 Application to 'real' radio sources 



p irts 



(if 



The cocoons of radio sources are elongated along the jet a: 
cause of the different external density the various source 
immersed in, this deviation from purely spherical shape ma) 
a different dynamical behaviour at the tip of the cocoon co npared 
to its base. In steady-state this adds only numerical factor 
solutions discussed above. However, it may prolong the trinsition 
to the coasting phase at the end of the lifetime of the jets 
vestigate this possibility, fully 3-dimensional simulations 
sour ces making the transition from active to coasting are 
(e.g. Reynolds, Heinz & Begelman, 2001). To date the e 
from these simulations is somewhat inconclusive. Neverth 
is unlikely that the geometrical shape of the radio source 
influences its evolutionary behaviour. 

So far we have implicitly assumed that the information that 
jets have ceased to supply energy to the cocoon propagates 
taneously through the whole cocoon. This is a good approx 
if the cocoon is small and the internal sound speed is 

(199"/| ) show that in the absence of 



veiy 



Kaiser & Alexander 



cant mixing of the cocoon material with the surrounding ga: 
the cocoon boundary the internal sound speed is indeed 1 
cases where mixing is important, the sound travel time al 
cocoon may become comparable to or even exceed the dyi lamical 
timescale. Most of the cocoon will then continue for a sig lificant 
time to expand as if it was still ii^^he_active_^hase_(e^ Kaiser, 



tu dly 



Schoenmakers & Rottgering, 2000) 

In any scenario the pressure in the cocoon will even 
come comparable to the thermal pressure of the surrounding 
this point, the expansion will slow down and the bow shod 
ens. Once pressure equilibrium is reached, the expansion s 
together and the pressure does not change anymore. When 
this happens depends on the temperature of the external gas 
the transition of the cocoon from expansion to pressure equifibrium 
will be a gradual process with those regions of the cocoon 
the centre of t he density distri bution coming into pressure 
rium first (e.g. Baldwin, 1982). Even in this equilibrium pi 



cocoon will not be static. As it is underdense with respec 
surrounding gas it will start rising due to buoyancy (e.g. 
Northover, 1973 



Churazov et al., 2001 



Brtiggen & Kaisei 



and the pressure inside the cocoon is decreasing. Howe 'er, the 
timescale for this rise usually exceeds the lifetime of the rel, .tivistic 
electrons. The equilibrium phase was also considered by |Komis- 
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sarov &. Gubanov (1994|) 
case of a radio source 
its surroundings at the 
also assume that the 
sure within the cocoon i 
In the following s 
or equilibrium radio sot 
cases. In each case the 
time /s < t, where / is 



For simplicity we study in this paper the 
w|hich comes into pressure equilibrium with 
ame time as the jets switch off. We will 
tralnsition is instantaneous and that the pres- 
5 uniform. 

udy of the radio emission from coasting 
rces we will concentrate on three limiting 
ets stop to supply energy to the cocoon at 
th 5 time of observation. 



Model A: The intern 
continues to evolve as 
i.e. i?c°=f-'''*^"'^' and p 

Model B: The intern 
immediately changes fn 
i.e. for Yc = Ys = 5/3 an 

Model C: The source 
roundings from the time 
between and t. 



il sound speed is slow and so the cocoon 
iflthe jets were still supplying it with energy, 
f(-4-P)/(5-P) even for f>/,. 
il sound speed is fast and so the cocoon 
■ )m the active phase to the coasting phase, 
1 / > /s we use equations (^). 
is in pressure equilibrium with its sur- 
fs. In this case R^^ and p are both constants 



3 RADIO EMISSIOIK FROM AN INACTIVE RADIO 
SOURCE 



The relativistic electron: 
chrotron emission are ac 
five jets and injected intc 
energy due to adiabatic 
ation and inverse Comp ti 
processes are cumulativ 
the cocoon later may 
may have faded already 



Alexander (1997) for 
of the energy distributic 
given time The total 
ing up the contribution; 
between f; = and f; = 
In the fo llowing we wi f 
tailed in 



Kaiser (2000) 



sources in the coasting 

The model employ 
having power-law depei d 
vious sections that this 
coasting phase. We find 
Model A by simply setltn: 
fi to /s < f at which the 
B and C w e have to sli; ;htlv 



et al. (1997). The relati 
time f; during a short tirr|e 
To determine the radio 
know the energy distribiition 
ergy density of the mag: 
determined by the evolu 
of time. Defining 6Vs = 
at t > h 



t \ Tc(7+3Tc-2P 



6(Tc + 



8V = 5Vs 



where we have assumec. 
adiabatic expansion 
sum of all volume elem 
of the cocoon. The time 



responsible for the emission of radio syn- 
elerated by the shocks at the end of the ac- 
the cocoon. They subsequently loose their 
;xpansion of the cocoon, synchrotron radi- 
on scattering of CMB photons. All these 
; and so relativistic electrons injected into 
ill radiate while material injected earlier 
The model of Kaiser, Dennett-Thorpe & 



a( tive 



sources allows to track the evolution 
n of the relativistic electrons injected at a 
radio spectrum is then found by integrat- 
of all electrons injected into the cocoon 
, where t is the current age of the source, 
use this model with the modifications de- 
ilso for the calculation of the emission of 
ase. 

;d relies on all relevant physical quantities 
lencies on time. We have seen in the pre- 
s, at least approximately, the case for the 
the total radio emission of the cocoon for 
g the upper limit of the integration over 
jets switched off. In addition, for Models 
modify the basic approach of K aiser 



'istic electrons injected into the cocoon at 
interval 8t[ are contained in a volume SV. 
ynchrotron emission from 6V, we need to 
of the relativistic electrons and the en- 
: letic field inside 8V. These are completely 
ion of 8V and the pressure, p, as functions 
5V (is) and ps= p{ts), yields for Model B 



and p = ps 



(8) 



that each volume element 8V undergoes 
1 5V o= i?^. This follows naturally as the 
;nts 8V must be equal to the total volume 
dependencies are then given by equations 
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Table 1. Model parameters for the fiducial model (see text). 



fixed parameters 


Rl 

s 

Ymax 
P 
Z 


2 

2.5 
10* 
1.5 
0.1 






1 1 years 


10^ 


c 
o 


free parameters 


fs / years 


10^ 






Ps 1 ergs cm"^ 


3 X 10"" 




deiived parameters 


2j / ergs s" ' 

P / -IS 
Poa(;/gcm 


1046 

10' 





Obviously, for Model C we have 8V = 8Vs and p = ps for all 
times f > /s- With this it is now straightforward to compute the syn- 
chrotron emission of each 8V and numerically in tegrate over the 



injection time t\ as desc ribed in Kaiser et al. (1997) 



Kaiser et al. (199 /) assume efficient pitch-angle scattering for 
the relativistic electrons moving in the magnetic field. Th is is usu- 
ally referred to as JP-type models (Faffe & Perola, 1973). For in- 



efficient scattering, the electrons with large pitch-angles, 6, loose 
their energy faster than those with small 6 (KP-type models: K ar- 
dashev, 1962; pacholczyk, 1970| ). It is straightforward to modify 
the present model to account for negligible pitch-angle scattering 
by assuming a uniform pitch-angle distribution and the introduc- 
tion of one further integral over 6. We calculate the radio emission 
from all models with and without pitch-angle scattering. 

The model for the radio emission depends on th e following 
parameters (see also 



Kaiser et al., 1997 



Kaiser, 2000): t the cur- 



rent age of the source, the time at which the jets switched off, ps, 
the pressure in the cocoon at ts,RT, the aspect ratio of the cocoon, s, 
the exponent of the initial power-law energy distribution of the rel- 
ativistic electrons, Ymax, the high energy cut-off of this distribution, 
P, the exponent of the density distribution of the external m edium 
and, z, the cosmological redshift of the source. Kaiser (200C ) shows 
that most of these parameters are degenerate in the sense that their 
values do not determine the shape of the radio spectrum in a unique 
way. We therefore concentrate on the effects of changing ps, t and 
fs and keep all other parameters fixed at their fiducial values given 
in table 111 We also assume that, at least at the time of injection 
into the cocoon, the magnetic field is in equipartition with the rel- 
ativistic particles in each 8V. We neglect any contribution to the 
cocoon pressure from non-radiating particles, i.e. thermal electrons 
or protons. With these assumptions and free parameters a model 
spectrum is completely determined except for a multiplicative nor- 
malisation constant, Af. For the fiducial model we set Af = 1. From 
Ps, ts and Af it is possible to infer the energy transport rate of the 
jets, Qj, and Poflg, which descr ibes the dens ity distribution of the 
gas surrounding the source (see Kaiser, 2000, for details). 

For the fiducial model as described in table |l|, f = f.s and so 
the jets are still active. The curves labeled 'f in Figure ^ show 
the spectrum for this model with (JP) and without (KP) pitch-angle 
scattering. With an age of 10^ years, the fiducial model predicts 
a fairly aged synchrotron spectrum. However, the cumulative en- 
ergy losses of the relativistic electrons due to synchrotron radiation 
and inverse Compton scattering are more severe at high frequency 
for the JP-type model. The other model spectra in Figure |^ are for 
the same model parameters as the fiducial model but with an age 
f = 2 X 10^ years which implies a coasting phase of 10^ years after 
fs. The spectral shape is very similar among all JP-type and KP-type 




100 1000 10000 

Observing frequency in MHz 

Figure 3. Model results for various coasting radio sources. Solid lines: With 
efficient pitch-angle scattering, dashed lines: Without pitch-angle scatter- 
ing. Labels refer to the various models defined in the text, 'f ' is the fiducial 
model and curves labeled 'A', 'B' and 'C are for the same source but with 
r = 2 X 10* years. 



models so that they appear to be the same spectrum simply shifted 
along the two coordinate axes. This is caused by the, for the chosen 
parameters, rather weak energy losses of the relativistic electrons 
after fs. The break in the energy distribution of these electrons is 
set mainly by the losses before fs and so occurs at roughly the same 
energy in all models. The break in the spectrum resulting from this 
is then simply shifted in frequency/flux space by the strength of 
the magnetic field. Since the energy density of the magnetic field is 
proportional to the pressure in the cocoon, it follows from Section 
that models A and C should have the lowest and highest break 
frequency respectively. This is confirmed by Figure bI Note here 
that for our choice of Yc and P the behaviour of p as a function of 
time for models A (p oc f-i'/?) and B (p oc f-i*'/'') is very similar, 
even for cases with larger energy losses of the relativistic electrons. 
In fact, usually < P < 2 and as long as Yc = 5/3, Models A and B 
are always similar. As expected, the spectral break occurs at higher 
frequencies for the KP-type models compared to the JP-type mod- 
els. 

Figure ^ shows a comparison of the effects of varying the free 
parameters f, fs and ps. For this comparison we use the fiducial 
Model A with effective pitch-angle scattering and f = 2 x 10^ years. 
Since the energy densities of the magnetic field and that of the rel- 
ativistic electrons is proportional to the pressure in the cocoon, an 
increased ps results in a higher radio flux at low frequencies. How- 
ever, the energy losses of the particles due to synchrotron emission 
are also increased moving the break in the spectrum to lower fre- 
quencies. Keeping the ratio fs/f constant but decreasing the source 
age to f = 10^ years, results in a brighter source at low frequencies 
with the spectral break moving to higher frequencies. At low fre- 
quencies the younger age implies a higher cocoon pressure which 
leads to an increased radio flux. As the radiative energy losses of the 
relativistic electrons are cumulative, the younger source age simply 
translates to smaller total losses and therefore a higher break fre- 
quency. 

Finally we note here that the injection of the relativistic elec- 
trons by the jets into the cocoon over a range of time leads to a 
flatte r spectrum at high frequencies than for a single burst injection 
(e.g. Leahy, 1991). 'Mild' breaks observed in the spectra of ra- 



dio sources are often interpreted, when using single burst injection 
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Figure 4. Model results for varying model parameters. Solid curve: Model 
A for fiducial parameters but with ; = 2 x 10** years and efficient pitch-angle 
scattering (same as solid line labeled 'A' in Figure Short-dashed curve: 
Same model but with = 1.5 x 10"'^. Long-dashed curve: Same model 
but with f = 2 X 10^ years and = 10^ years. 



models (as assumed in Komissarov & Gubanov, 19941, as result 



ing from the suppression of pitch-angle scattering. For a range in 
the injection time, a similar flattening may be achieved even for 
efficient pitch-angle scattering. 



4 COMPARISON WITH THE SPECTRA OF RADIO 
RELICS 

In this section we will use the models developed above to fit the ob- 
served spectra of radio relic sources. In order to compare our results 
with those of Komissarov & Gubanov (1994), we use their data 



points and error estimates for the same relic sources, i.e. Mol 0100- 
22 (orMRC 0100-221), Cul 0038-096, 4C 38.39, 4C 63.10, 3C 318 
and 3C 464. The names of the clu sters these sources are found in i s 
given in table 0. As discussed in Komissarov & Gubanov (1994), 



there is a contr adiction between the flux measurements of S lee & 
Reynolds (1984) and |Slee, Perley & Siegman (1989|) at 1.46 5 GHz 
for Cul 0038-096. We follow [Komissarov & Gubanov (1994| ) in us- 
ing both flux measurements separately, assigni ng indices 'a' and 'b ' 
to the different spectra. However, we note that Bagchi et al. (1998) 
and t he relic flux density meas ure from the NVSS are in agreement 
with Slee & Reynolds (1984). We also fit our models to the ra- 



dio spectrum of the relic source 1253-1-275. The relic is somewhat 
different to the six others in that it is rather extended and shows 
practically no break in its radio spectrum. Because of these differ- 
ences, 1253-1-275 may have had a completely different origin and/or 
history compared to the other relics. However, we will show below 
that the radio spectrum of this source is also well fitted by our mod- 
els. The ra dio fluxes for 1253-1-275 at seven radio frequencies can 
be found in Giovannini et al. (1991 ) 



We fit the spectra adopting the fixed model parameters of the 
fiducial model and using t, and ps as free parameters. The good- 
ness of the fit is determined by the standard x^-technique. The nor- 
malisation of the spectrum, Af, is chosen such that for each set of 
free parameters the x^-value is minimised. The minimum of as 
a function of the free parameters is found using a d ownhill simplex 
method in three dimension (e.g. Press et al., 19921. To avoid solu- 
tions with excessively young or old sources, we restricted t to the 
interval from 5 x 10^ years to 5 x 10** years. 



Figure 5. Model fits to the observed spectrum of Mol 0100-221. Subscripts 
'JP' denote models with efficient pitch-angle scattering, subscripts 'KP' 
are assigned to models without pitch-angle scattering. Flux measurements 
are shown with their errors. Dotted line: Model Ajp. Short-dashed line: 
Model Bjp. Dash-dotted line: Model Cjp. Dash-triple dotted line: Model 
Akp- Long-dashed line: Model Bjcp. Solid line: Model Ckp- 
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Fi gure 6. Same as Figure H but for Cul 0038-096 using flux measurement 
ot isiee& Reynolds (1984 at 1.465 GHz. 
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Fi gure 7. Same as Figure |^ but for Cul 0038-096 using flux measurement 
of fice etal. (198s| ) at 1.465 GHz. 
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Table 2. Best fitting model parameters. The Abell number of the host cluster of the respective radio rehc is given below the source name. Subscripts ' 
denote models with efficient pitch-angle scattering, subscripts 'KP' are assigned to models without pitch-angle scattering. 
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Figure 8. Same as Figure H but for 4C 38.39. 
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Figure 10. Same as Figure | but for 3C 3 1 8. 1 . 




100 1000 
Observing frequency in MHz 



Figure 9. Same as Figure H but for 4C 63. 10. 
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Figure 11. Same as Figure ^ but for 3C 464. 



The best-fitting model parameters are summarised in table ^ 
and the resulting spectra are plotted in figures I through In} In 
practically all cases the best-fitting spectra from the various models 
are almost indistinguishable from each other. In the following we 
briefly discuss the results for each source individually. 

Mol 0100-221: The spectrum is well fitted by all models (re- 
duced between 0.99 and 1.03). The JP and KP-type models give 
somewhat different results for the free model parameters. How- 
ever, within each family of models they are roughly consistent with 
each other. This indicates that the shape of the spectrum is mainly 
determined during the active phase of the source. The exact be- 
haviour of the cocoon in the subsequent coasting phase has little 
influence. The JP-type models prefer older ages, higher jet pow- 
ers and less dense environments compared to the KP-type models. 
KP-type models are marginally preferred over JP-type models. 

Cul 0038-096: Both variants of the source spectrum can be 
successfully fitted by the models. However, the goodness of fit and 
the values of the free parameters vary significantly for each variant. 
The best fit achieved has a reduced -value of only 0.3, but pre- 
dicts a very low density of the external gas. Three of the bets-fitting 
models (Akp for variant a, Akp and Bkp for variant b) have ages 
very close to the maximum age, f = 5 x 10^ years, allowed for the 
fit. At least for models A and B the JP-type and KP-type models 
give similar results, but models C are quite different from these. 



The spectrum of this source can be fitted reasonably well with all 
the models discussed here. However, a clear trend for the model 
parameters is not apparent. 

4C 38.39: For this source models with high external densities 
are clearly preferred. However, the spread in the value of poag is 
considerable. There is also no agreement between the models with 
respect to the source age. From figure ^ it seems that the rather flat 
spectral break at high frequencies causes problems for all models. 
Again there is no discernible trend in the best-fitting parameters. 

4C 63.10: The spectrum of this source is fitted about equally 
well by all models. In terms of the jet power, Qj, and the external 

density parameter, Poag, JP and KP-type models agree well with 
each other for models A and B, but there is a large range for the 
predicted source age. Models C require a somewhat lower density 
and they predict very low source ages. 

3C 318.1: The model fits to the spectrum are very good with a 
reduced of around 0.12. Values for Qj and PoAq are similar for 
all models but again for this source there is a large spread for the 
source age. 

3C 464: Except for model Ajp, all models predict rather low 
source ages. The best fits are achieved by both models C with rather 
high external densities. The predicted values of Qj are fairly uni- 
form but the external density vary strongly between models. 
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Figure 12. Same as Figure |but for 1253+275. 

1253+275: The model parameters agree with each other for 
KP-type and JP-type models. However, the parameters are incon- 
sistent between the two groups. The KP-type models prefer a young 
source age with virtually no influence of the phase in the source 
evolution when the jets are switched off. The predicted densities 
and pressures are very high. In the absence of radiative energy 
losses,our choice of p = 2.5 for the power law exponent of the en- 
ergy distribution of the relativistic electrons at injection time, im- 
plies a spectral index of about —0.75. The observed spectrum of 
1253-1-275 shows a spectral index of around —1.1 with no spec- 
tral break. The in general flatter spectra resulting from the KP- 
type models do not fit this steep spectrum well. Therefore, the 
model requires high external densities causing high pressure inside 
the lobes which subsequently show substantial energy losses and 
steeper spectra. The predictions of the JP-type model for an old 
source with moderate pressure values are more realistic. The low 
jet powers predicted by all models are caused by the substantially 
lower radio flux of this relic compared to the other sources. 



5 CONCLUSIONS 

We have developed a model for the evolution of the lobes of radio 
galaxies and radio-loud quasars of type FRII after the jets of the 
central AGN have stopped supplying them with energy. In the case 
that the lobes are still overpressured with respect to their surround- 
ings, we show that their volume is proportional to 

6(Tc + l) 

;forYc =4/3 andYs =5/3. 



V oc 



(9) 



The proportionality given above for yc = Ys is, strictly speaking, not 
a solution of the equations governing the hydrodynamics. However, 
we have shown that is provides an excellent fit for the numerically 
found correct solution. The numerical calculations also show that 
the transition from the active phase with jets to the coasting phase 
without jets is fast. 

We identify the lobes in the coasting phase with radio relic 
sources. The calculated model spectra depend only on a small num- 
ber of source and environment parameters: The total age of the 
source, t, the age of the source when the jets were switched off, 
ts, and the pressure inside the lobes at that time, ps- We consider 
three limiting cases for the lobe evolution: The lobes continue to 
expand as if the jets were still supplying energy (Model A), the 



lobes continue to expand but with the modified dynamics found 
for the coasting phase and summarised above (Model B) and the 
lobes come into pressure equilibrium at the time the jets switch 
off and stop expanding (Model C). We also investigate the effects 
of efficient pitch-angle scattering (JP-type models) and its absence 
(KP-type models). 

All our models can provide satisfactory to excellent fits to the 
observed spectra of radio relics. Taking into account the evolution 
of the strength of the magnetic field inside the lobe during the in- 
jection of relativistic electrons flattens the spectra sufficiently to ex- 
plain the observations. Thu s we cannot rule out efficient pit ch-angle 
scattering as suggested by Komissarov & Gubanov (1994). Unfor- 



tunately, the good fits provided by all models also imply that the 
radio spectra of relic sources alone do not constrain any of the im- 
portant source or enviro nment parameter s. The situation is aggra- 
vated by the findings of [Slee et al. (2001 ) that the spectra can also 
be fitted using the alternative assumption of inhomogeneous mag- 
netic fields inside the relics. Furthermore, even the inactive lobes 
may be re-ignited by c ompression during cluster mergers ( EnBlin 
& Gopal-Krishna, 2 001). Note here that the model of E nBlin & 
Gopal-Krishna (2001) for radio relics does not require the radio 
plasma to originate in a powerful FRII-type radio galaxy. 

The morphologies of radio relics are very varied and only 
rarely resemble the well-defined lobes of active radio galaxies. This 
strongly suggests that fluid flows in the clusters significantly influ- 
ence the evolution of radio galaxy remnants. Any re-acceleration or 
energisation due to compression would make it even more difficult 
to determine the original conditions in the relic. A possible solution 
could be the combination of radio spectra with X-ray observations 
of sufficient spatial resolution to identify inverse Compton scattered 
photons originating in the relativistic plasma of the relics. The X- 
ray observations could then be used to estimate the electron density 
while the radio emission would constrain the strength of the mag- 
netic field. Despite the problems, it is interesting to note that our 
simple models can fit the radio spectra of all relic sources stud- 
ied here. This includes the clearly curved spectrum of Cul 0038- 
096, a relic with substantial substructure, as well as the spectrum 
of 1253-1-275 without any significant spectral break and a compara- 
tively smooth appearance. This would imply that, whatever the his- 
tory of the relic source, even a rather simplistic model can explain 
the observed spectra in a variety of different ways thus making it 
impossible to constrain the conditions within the relic sources. We 
cannot rule out a complex history or re-acceleration of relativistic 
electrons influencing the properties of the relic sources. However, 
in the light of our results, such more complicated models are not 
required to fit the available radio data. 

The degeneracy of model parameters can, however, also be 
turned to advantage in the case of statistical studies of the ra- 
dio source population. For example, we can use the known radio 
luminosity function of active FRII objects and make predictions 
for the cosmological distribution of their remnants. Because of 
the mentioned degeneracy, these predictions will be almost model- 
independent. Such statistical studies can be used to predict the de- 
tection rate of radio galaxy remnants in low frequency radio sur- 
veys with high surface brightness sensitivity (Cotter & Kaiser, in 
preparation). 

In fact, it is not yet clear whether radio relic sources are ac- 
tually the remains of the lobes of powerful radio galaxies. Deep 
observations of radio galaxies in which the jet flows have stopped 
comparatively recently are necessary to decide whether the con- 
ditio ns inside their lobe s point towards them evolving into relics 
(e.g. Kaiser et al., 200C ). Clearly, the radio relic sources continue 
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to challenge our understanding of the physics of radio galaxies and 
clusters. 
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